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SUMMARY 


A  continuum  mechanics  description,  of  the  phenomenon  of  stress  Induced 
microcracking  has  been  used  to  study  the  near  tip  stress  and  strain  fields  and 
the  size  and  shape  of  a  small  scale  damaged  zone  for  a  stationary  mode  I  crack 
in  an  elastic  body.\  The  material  model  is  characterized  by  a  microcracking 

criterion,  which  is  an  extension  and  simplification  of  the  generalized 

/ 

microcracking  criteripn  proposed  by  Fu  and  Evans  [10]  for  the  case  of  thermal 
stress-induced  micyocracking.  That  together  with  approximate  expressions 
relating  the  effective  composite  moduli  to  the  elastic  properties  of  the 
brittle  material,  via  the  microcrack  density  e  first  introduced  by  Budiansky 
and  O'Cdnnell  [1],  yield  a  self  consistent  approach  to  the  stress  induced 
micr^rracklng  phenomenon. 

>/ The  microcrack  density  was  found  to  characterize  three  regions  of 
interest.  In  the  outer  region  the  microcrack  density  is  zero  and  the  stress 
and  strain  fields  are  purely  those  for  linear  elastic  deformation.  This 
elastic  field  constrains  the  microcracking  deformation  which  in  combination 
with  material  weakening  due  to  microcracking  causes  stress  relaxation  in  a 
region  of  intermediate  microcracking.  Very  close  to  the  crack  tip  the 
microcrack  density  is  saturated  and  the  stress  field  becomes  again  singular 
but  with  a  lower  stress  intensity  than  would  prevail  in  the  absence  of 
microcracking.  In  the  case  where  very  rapid  microcracking  occurs  as  the 
strain  is  increased,  the  intermediate  microcracking  zone  is  still  present 
providing  continuity  of  the  strain  field  and  a  smooth  transition  of  the  stress 
field  from  the  purely  elastic  region  to  the  region  with  saturated  microcrack 

density. It  appears  that  the  existance  of  the  region  of  intermediate 

\  '  -N 

microcrack  density  is  essential  to  preserving  the  assumption  of  a  continuum 
composite  be  '  se  it  helps  to  avoid  any  strain  and  stress  discontinuity. 


INTRODUCTION 


The  self  consistent  aproach  to  the  stress-induced  microcracking  of 
brittle  materials,  as  already  mentioned.  Is  based  on  two  major  notions*  The 
notion  of  a  microcracking  criterion  and  that  of  effective  elastic  properties 
of  the  mlcrocracked  material*  The  microcracking  criterion  relates  the  induced 
microcrack  density  e  to  the  existing  stress  magnitude.  The  effective  elastic 
properties  of  the  mlcrocracked  material  are  then  functions  of  the  Initial 
elastic  properties  and  the  microcrack  density.  Budansky  and  O’Connell  [1], 
Introduced  the  microcrack  density  parameter  e  given  by  equation  (1) 


e 


(1) 


where  N-number  of  microcracks  per  unit  volume,  A  is  the  area  of  the  microcrack 
surface  and  P  Is  the  perimeter  of  the  microcrack.  Using  energy  balance 
considerations  and  fracture  mechanics  analysis  they  were  able  to  derive 
mathematical  expressions  relating  the  effective  elastic  properties  of  the 
actual  state  to  those  of  the  unmlcrocracked  elastic  state,  through  the 
microcracking  density  e.  We  will  make  use  of  the  work  of  Fu  and  Evans  [10]  to 
construct  a  microcracking  criterion  which  will  relate  increases  in  the 
microcracking  density  to  increases  In  stress  magnitudes*  Together  with  the 
self  consistent  results  of  Budlansky  and  O'Connell  [1],  the  microcracking 
criterion  can  be  used  to  develop  a  constitutive  law  for  the  inelastic  behavior 
of  a  microcracking  continuum. 

Our  Interest  Is  In  the  development  of  microcracks  near  the  tips  of  major 
cracks.  To  study  this  phenomenon  In  an  approximate  way,  we  use  the 
constitutive  law  for  the  microcracking  material  In  finite  element  calculations 
of  either  plane  strain  or  plane  stress  deformation  near  the  crack  tip.  These 
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results  are  relevant  to  the  question  of  mlcrocrack  toughening  In  brittle 
materials.  Before  we  proceed  to  describe  these  results  we  will  discuss 
briefly  some  previous  work  on  this  problem. 

Hoagland,  et  al.  [2],  proposed  a  simplified  method  for  estimating  the 
density  of  microcracks  In  the  vicinity  of  a  major  crack  tip.  They  assume  a 
set  of  randomly  oriented  lines,  which  are  the  traces  of  the  potential 
mlcrocrack  planes,  and  which  become  microcracks  when  the  stress  normal  to  the 
line  exceed  a  critical  value.  They  also  assume  that  the  dimension  of  the 
microcracks  Is  very  much  less  than  the  major  crack  size  of  the  model.  In 
determining  the  stress  Intensity,  they  used  the  singular  elastic  stress  field 
near  the  crack  tip.  Later  studies  done  by  Evans  [3],  suggest  that  the  uniform 
mlcrocrack  density  analysis  presented  by  Hoagland  et  al.  is  inaccurate  because 
there  is  an  increased  density  of  microcracks  at  the  higher  stress  levels, 

which  more  strongly  biases  the  mlcrocrack  density  towards  the  crack  tip. 

In  L980  Hoagland  and  Embury  [4],  suggested  a  numerical  procedure  for 

modeling  the  number  and  distribution  of  microcracks  around  a  crack  tip,  as  a 

function  of  the  applied  stress  intensity.  The  procedure  accounts 

approximately  for  microcrack-microcrack  and  microcrack-crack  Interactions. 
This  work  Is  actually  an  extension  of  their  original  concept  of  nucleating 
microcracks  based  on  a  uniform  fracture  stress  by  adding  new  stress  fields, 
whenever  a  new  mlcrocrack  is  nucleated.  Although  Hoagland  and  Embury  arrived 
at  some  useful  conclusions  concerning  the  stress  induced  microcracking  In  the 
vicinity  of  a  major  crack,  their  results  are  restricted  to  a  two-dimensional 
problem  with  through  thickness  microcracks  extended  along  the  third 
dimension.  In  contrast  to  that  the  theory  developed  by  Budlansky  and 
O'Connell  accounts  for  microcracks  with  finite  third  dimension  (i.e.. 


circular,  elliptical,  etc*)  and  allows  this  effect  to  be  Introduced  in  an 
approximate  way* 

A  more  thorough  investigation  of  the  mechanics  of  microcrack  toughening 
was  carried  out  by  Fu  and  Evans  [10]*  In  their  work  they  established  the 
basis  of  a  continuum  mechanics  description  of  the  microcrack  toughening 
process*  The  constitutive  behavior  of  the  microcrack  zone  is  based  on  the 
overall  strain  response  of  a  microcracking  medium,  subject  to  a  uniform  stress 
field.  A  generalized  microcracking  criterion  for  microcracking  at  facets 
subject  to  general  stress  is  proposed.  Stress  Induced  microcracking  is 
initiated  at  a  threshold  load.  Above  the  threshold  load,  the  stress/straln 
behavior  becomes  nonlinear.  During  unloading,  linearity  resumes  and 
hysteresis  is  present,  corresponding  to  dissipation  of  strain  energy. 

THE  MICROCRACKING  CRITERION  AND  CONSTITUTIVE  LAW 

Fu  and  Evans  [10]  studied  the  phenomenon  of  stress  induced  microcracking 
of  grain  boundary  facets.  These  facets  are  generally  subject  to  residual 
stresses  caused  by  thermal  expansion  anisotropy  of  the  neighboring  grains. 
When  a  high  tensile  stress  is  applied,  favorably  oriented  facets  crack.  The 
crack  is  confined  to  Chat  facet  and  it  takes  a  much  higher  stress  to  cause  it 
to  propagate  out  of  its  initial  grain  boundary  interface.  As  a  consequence, 
the  process  of  microcracking  is  at  least  initially  stable,  with  a  steadily 
increasing  stress  causing  the  nucleation  of  more  cracks  rather  than  the 
propagation  of  existing  ones. 

From  experiments  and  fracture  mechanics  analysis,  Fu  and  Evans  concluded 
Chat  a  particular  facet  will  crack  if  it  is  larger  than  a  critical  size  which 
depends  on  the  resolved  applied  stress  and  the  residual  stress  on  the 
interface.  In  this  notion  is  implicit  the  idea  that  facets  can  microcrack  due 


Co  the  residual  stress  alone  if  the  facet  exceeds  a  critical  size*  Fu  and 

Evans  were  able  to  phrase  their  results  in  terms  of  a  parameter  e  Introduced 

by  Budiansky  and  O'Connell  [1]  to  describe  microcrack  density*  When  there  are 

N  circular  microcracks  per  unit  volume  and  the  microcracks  are  of  uniform 

3 

size  with  radius  a  then  e  ■  Na  .  Fu  and  Evans  studied  microcracking  due  to 
macros coplcally  uniform  biaxial  stress  fields  (with  principal  stresses 
0j>02,  k**02/oj).  From  their  analysis  and  observations,  they  proposed 

e  -  k(o)(o1  -  crc)  for  >  ac  (2) 

where  X  is  a  parameter  that  depends  on  stress  state  and  material  properties 
but  independent  of  stress  magnitude.  For  example  X  for  k-1  is  twice  the  value 
of  X  for  k-0.  The  parameter  e  represents  the  density  of  microcracks 
introduced  by  the  applied  stress  that  is  in  excess  of  any  initial  density. 
Thus  the  microcrack  density  e  Increases  linearly  with  stress  and  does  so  up  to 
a  saturation  level,  after  which  it  remains  constant.  e  is  unchanged  by  a 
reduction  of  stress. 

We  require  a  criterion  for  microcracking  which  can  be  used  in  more 
general  states  of  multlaxlal  stress.  For  this  purpose  we  will  modify  equation 
(2)  to  become. 

e  -  X(oR  -  oc)  if  oR  >  oc  (3) 

where  oR  -  is  an  effective  stress,  <Jq  is  the  critical  stress  for 

microcracking  initiation  and  X  is  a  material  constant  independent  of  stress 
state.  Of  course  equation  (3)  only  applies  when  the  largest  principal  stress 
is  tensile  and  Increases  monotonlcally.  The  criterion  (3)  is  somewhat 
different  from  that  proposed  by  Fu  and  Evans  [10],  If  we  take  the  uniaxial 
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stress  case  with  ■  o  as  a  datum  then  both  criteria  give  e  “  Xu(o  ”  °g) 
where  X^  is  the  value  of  X  from  Fu  and  Evans  theory  for  uniaxial  tension.  For 
equal  biaxial  tension  (o^  ■  m  o),  Fu  and  Evans'  criterion  gives 
e  *»  2X^(0  -  oc>  whereas  the  one  we  have  proposed  leads  to  e  ■  Xu(/2o  -  o^,)  . 
Thus  our  criterion  shows  a  reduced  rate  of  microcracking  In  biaxial  tension 
compared  co  the  rate  suggested  by  Fu  and  Evans  but  microcracking  starts  at  a 
lower  stress  according  to  our  criterion.  It  can  be  argued  that  the  microcrack 
density  in  biaxial  tension  would  not  be  double  that  for  uniaxial  stress,  since 
a  given  facet  cannot  be  cracked  twice.  That  is,  unless  the  families  of  facets 
cracked  by  each  applied  stress  component  are  mutually  exclusive.  However  we 
are  interested  mainly  in  a  simple  constitutive  law  with  appropriate 
characteristics  for  microcracking  to  illustrate  what  will  occur  near  major 
crack  tips.  We  will  develop  more  exact  microcracking  laws  in  future  work. 


We  summarize  our  microcracking  criterion  as  follows. 
For  monotonlcally  Increasing  stress,  if 


<  Oq  then  €  -  0 


Material  remains  unmlcrocracked 


°C  ^  aR  <  °M  then  e  *  ”  °C)  The  raicr°crack  density  is 

increasing  linearly  with  a 


Chen  e  -  X(0jj  ~  °c)  *  ^  The  microcrack  density  is  saturated. 


Furthermore  e  cannot  decrease.  The  stress/strain  relation  given  for  a 
macroscopic  element  of  material  is  linear  for  non-mi crocracked  regions  (e**0). 


becomes  nonlinear  for  regions  with  increasing  microcrack  density  0<e<e^  and 
again  becomes  linear  for  fully  microcracked  regions  (e  *  e^).  During 
unloading  linearity  resumes  and  when  material  unloads  to  zero  stress  no 
permanent  strains  are  present  as  shown  in  fig.  (1). 


Nv*--3 


-  Jk  »  •  J  *•  kj»*  *  .  *  L  •  o'  >  ■ -  r  >> 


IC  should  be  noted  that  the  alcrocrack  distribution  produced  by  a 


deviatoric  applied  stress  coupled  with  the  residual  stress  will  be 
anisotropic.  However,  Fu  and  Evans  used  the  isotropic  theory  of  Budiansky  and 
O'Connell  [1]  to  determine  macroscopic  moduli  for  the  mlcrocracked  material, 
that  is  they  used  e  as  if  it  represented  a  random  distribution  of  microcracks 
rather  than  an  oriented  one.  We  will  follow  Fu  and  Evans  in  this  regard  to 
obtain  an  approximate  continuum  theory  for  mi croc racking. 

We  have  found  that  we  can  approximate  Budiansky  and  O'Connell's  results 


as  follows. 


E  v  .  16  1 

mi  _  m  1  —  ■  p  m  _ 

E  v  9  f 


where  E,  v  and  E,  v  are  Young's  moduli  and  Poisson's  ratio  for  the  material 
before  and  after  stress  Induced  microcracking  respectively.  Note  that  E  and 
v  are  defined  at  constant  £.  Thus  the  constitutive  equation  for  the 
microcracking  solid  becomes 


E+v  „  v  , 

E  °ij  t  *kk  6ij 


where  g,  is  the  macroscopic  strain,  a  is  the  macroscopic  stress,  is  the 

kronecker  delta  and  rules  for  f  are  given  by  equations  (3),  (4)  and  (5). 

The  constitutive  equation  (6)  is  only  valid  for  values  of  the 

microcracking  parameter  f  within  the  interval  l<f<®,  fig  (3).  This  together 

with  equation  (5)  defines  the  limiting  values  of  the  microcracking  density  £, 

9 

0  <  £  <  (See  fig  (3)).  The  upper  bound  of  £  with  the  aid  of  equation  (3) 
yield  a  relation  between  the  critical  stresses  a a M  and  the  parameter  X 
i.e.,  °M  "  ac  <  JZ  *  X  (See  fi« 


FORMULATION  OF  CRACK  TIP  PROBLEM 


The  high  stresses  near  the  tip  of  a  long  crack  will  cause  microcracking 
in  the  near  tip  region  in  the  material  of  Interest  to  us,  as  shown  in  fig 
(5).  We  will  restrict  our  attention  to  the  situation  where  the  zone  of 
microcrack  damage  is  very  small  compared  to  the  body  containing  the  crack.  We 
will  refer  to  this  situation  as  small  scale  microcracking*  In  this  case,  the 
microcrack  zone  will  lie  within  a  region  of  uncracked  material,  some  distance 
outside  the  zone,  the  stresses  will  be  almost  the  same  as  when  the  material 
does  not  microcrack.  When  small  scale  microcracking  prevails,  these  stresses 
will  be  the  singular  linear  elastic  crack  tip  stresses. 

It  follows  that  the  plane  small  scale  raicrocracking  problem  can  be  solved 
by  considering  a  plane  region  around  the  crack  tip  to  which  are  applied 
boundary  tractions  given  by  the  linear  elastic  stress  field,  i.e. 


no-  - 

J  ij  /HF 


=  °i  W9) 


where  T  are  the  boundary  tractions,  ri  is  the  outward  unit  normal  to  the 
region  boundary,  Kj  is  the  mode  I  (tensile  opening  mode)  stress  intensity 
factor  and  (r,6)  are  polar  coordinates  originating  at  the  crack  tip  as  shown 
in  fig  (6).  The  function  £  determines  the  angular  distribution  of  stress  and 

A* 

its  form  can  be  found  in  the  article  by  Rice  [23].  The  crack  surfaces  are 
traction  free.  To  ensure  small  scale  microcracking,  Kj  must  be  limited  to  a 
level  low  enough  that  the  zone  size  is  small  compared  to  the  region  for  which 
the  analysis  is  being  carried  out 

The  governing  equations  of  equilibrium  and  compatibility  are  enforced 


through  the  principle  of  virtual  work 
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JA  6etj  dA  -  Jg  Tt  6ut  ds  (8) 

In  the  absence  of  body  forces,  where  A  Is  the  plane  area  being  analyzed,  Sj  is 
the  perimeter  where  tractions  are  prescribed,  ^  are  the  displacements  and  the 
symbol  6  indicates  a  virtual  variation  of  the  quantity  following  it  and  the 
variation  disappears  on  S  -  Sj. 

The  constitutive  law  for  the  material  has  been  described  in  the  previous 
section.  However,  we  will  find  it  useful  to  state  the  following  form. 


Iij  =  Ci jkl  ekl  *  f+v  eij  +  Cf-2v)(f+v)  *kk  5ij 


(9) 


where  f 


and  e  =  e(o)  as  in  equation  (3) 


For  plane  problems 


*«P 


f+v 

E 


J  „  -  v*o  6  a] 
op  YY  <*P 


«,  P,  Y  *  It  2 


(10a) 


<xp 


E 

f+v 


v* 

<l-2v*)  eYY 


(10P) 


V  .  v 

where  v*  ■  -j  for  plane  strain  and  v*  *  for  plane  stress*  The  governing 
equations,  together  with  the  appropriate  boundary  conditions  corresponding  to 
a  specified  geometry,  define  a  boundary  value  problem  whose  solution  we  will 
obtain  numerically. 


THE  FINITE  ELEMENT  EQUATIONS 

The  finite  element  equations  can  be  derived  from  the  principle  of  virtual 
work  given  by  equation  (8).  Finite  element  interpolations  together  with 


equation  (8)  and  the  constitutive  law  (9)  and  (10.8),  give  rise  to  the 
nonlinear  finite  element  equations 

[k(un)]{un}  -  (F)  (11) 

where  [k]  is  the  stiffness  matrix,  (un)  the  nodal  displacement  vector  and  (F) 
the  force  vector.  The  notation  [k(un)]  Indicates  the  dependence  of  the 
stiffness  on  the  nodal  displacements  due  to  the  nonlinear  constitutive  law. 

In  the  finite  element  analysis,  a  4  noded  quadrilateral  isoparametric 
element  with  4  stations  for  the  integration  of  stiffness  was  used.  The  finite 
element  equations  (11)  were  solved  using  an  interative  method  updating  the 
stiffness  matrix  in  every  iteration.  The  microcracking  parameter  f  was  found 
to  be  the  root  of  a  sixth  order  polynomial  in  f  which  satisfies  the  conditions 
Imposed  on  f  earlier  in  this  work  i.e.  l<f<°°.  The  coefficients  of  the 

polynomial  are  strain  dependent  and  had  to  be  recomputed  in  every  iteration 
[see  Appendix  II]. 

To  ensure  small  scale  microcracking,  the  outer  radius  of  the  mesh  was 
chosen  to  be  15  to  20  times  the  outer  radius  of  a  damaged  zone,  estimated  from 
the  stresses  in  the  absense  of  microcracking. 

RESULTS 

Stress-strain  fields 

Figs  (7)  to  (8)  display  the  stresses  and  strains  ahead  of  the  crack  tip 
for  various  choices  of  ^  and  e^.  \  is  the  rate  of  increase  of  the  raicrocrack 

density  with  respect  to  the  effective  stress  and  eM  =  \(o^-oc)  is  the 

saturation  value  of  the  microcrack  density. 

In  the  unmicroc racked  region,  both  stresses  and  strains  agree  quite 

closely  with  the  solution  obtained  for  an  unmicrocracked  material.  The 

elastic  strain  and  stress  fields  constrain  the  deformation  due  to 


microcracking  in  the  damaged  zone.  The  above  constraint  together  with  the 
weakening  of  the  material  in  the  damaged  zone  due  to  microcracking,  causes 
stress  relaxation  in  the  region  of  intermediate  microcracking.  Very  close  to 
the  crack  tip  the  microcracking  density  is  saturated  and  the  stress  and  strain 
fields  become  singular  again.  We  observe  that  the  amount  of  stress  relaxation 
increases  with  X  and  it  reaches  its  maximum  level  for  \=*  ®.  At  this  point  we 
notice  the  importance  of  the  presence  of  the  zone  of  intermediate 
microcracking,  for  its  presence  ensures  a  smooth  transition  for  both  the 
strain  and  stress  fields  from  the  outer  to  the  inner  asymptotic  fields.  The 
absence  of  the  above  region  would  lead  to  infinite  strain  and  stress  gradients 
along  the  boundary  of  the  fully  microcracked  material  and  the  unmicrocrack 
material  [see  Appendix  I].  We  also  notice  that  in  the  stress  relaxation 
region,  the  effective  stress  which  in  our  case  is  the  equivalent  stress 
oR  *  / j  o^j  is  consistent  with  the  microcracking  law  eq  (3),  fig  (2)  that  we 
used.  In  the  particular  case  of  X=  ®,  oR=oc*1.0  as  we  expected.  In  this 
situation,  the  partially  microcracked  zone  is  like  a  perfectly  plastic  zone. 
The  zone  with  saturated  microcrack  density  Is  dominated  by  microcracking 
deformation  with  larger  strains  and  stresses  at  lower  levels  than  would 
prevail  ’  e  absence  of  microcracking. 


The  shape  and  size  of  the  transformed  zone 

A  typical  microcrack  zone  is  shown  In  fig.  (5).  Three  characteristic 
quantities  describing  the  microcrack  zones  that  we  studied  were  consistently 
computed.  xm  is  the  distance  along  the  x  axis  at  which  the  microcrack  zone 
boundary  crosses  the  x  axis.  h^  denotes  the  y  coordinate  of  the  farthest 
point  along  the  y  direction  to  experience  microcracking  and  9m  is  the  polar 
angle  corresponding  to  the  same  point.  The  boundary  of  the  inner  zone,  where 
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APPENDIX  II 


The  nicrocracking  law  given  by  equation  (3)  earlier  In  this  paper  allow 
us  to  determine  the  microcrack  density  e  for  a  given  stress  state*  It  Is  more 
convenient  to  determine  e  or  f  ■  1/(1  -  e)  for  a  given  strain  state. 


from  the  constitutive  law  [e^  +  eRR  ^ j  1 
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then  we  have  the  following  cases  in  computing  f. 
1)  aR(f)  <  ac  ->  e  *  0  ->  f  ■ 
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substituting  aR  from  equation  (2)  in  eqn  (5)  we  have 
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which  is  a  sixth  order  polynomial  in  f  of  the  form 


a^f^  +  a^f^  +  a^f^  +  a^f^  +  f^  +  a^f  +  Sq  *  0  (7) 


whose  coefficients  are  given  below. 
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stresa  discontinuity 


we  allow  axx  -  a%%  and  require  eyy  - 
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we  allow  e^  -  eyy  and  require  a ^  -  oxx 
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stress  and  strains  are  continuous,  for  mlcrocracked  materials  f>l  and 
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APPENDIX  I 


We  consider  an  element  consisting  of  two  dissimilar  material 


E,  v  are  the  Young  modulus  and  Poisson's  ratio  for  an  unraicrocracked  material* 
E,  v  are  the  corresponding  moduli  for  a  microcracked  material* 

We  require  that  the  material  remains  continuous  after  deformation.  Thus  from 
equilibrium  and  compatibility  we  get  the  interface  conditions 
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interface  conditions. 


The  constitutive  relations  for  plane  elasticity  case  are 
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where  e  is  the  microcrack  density. 


The  finite  element  analysis  will  substantially  contribute  to  a  better 
understanding  of  the  mechanics  of  stress  induced  microcrack  toughening  in 
brittle  materials  As  we  saw  earlier  in  our  analysis,  a  zone  with  saturated 
microcrack  density  always  exists.  We  may  be  able  to  treat  this  zone  as  being 
composed  of  homogeneous  elastic  material  and  use  a  Griffith  criterion  for 
propagation  of  the  major  crack  into  this  zone.  McMeeking,  in  unpublished 
work,  has  shown  that  the  material  obeying  the  constitutive  law  of  eq.  (3-6)  is 
nonlinear  hyperelastic.  It  follows  that  the  J-integral  is  path  independent 
throughout  the  material  in  both  unmicrocracked  and  mlcrocracked  regions*  As  a 
consequence,  the  energy  release  rate  for  the  microcracking  material  will  have 
the  same  value  as  in  the  nonmicrocracking  material  at  the  same  applied  loads. 
Thus  if  the  critical  value  for  propagation  is  unchanged,  then  no  microcrack 
toughening  can  occur.  However,  the  critical  energy  release  rate  for  fully 
mlcrocracked  material  may  differ  from  that  for  the  unmicrocracked  case.  A 
simple  model  of  the  crack  growing  from  microcrack  facet  to  microcrack  facet 
suggests  that  the  average  critical  energy  release  rate  would  be  less  and  so 
microcrack  embrittlement  would  occur.  However,  aspects  of  R-curve  type 
behavior  and  more  realistic  microcracking  laws  may  hold  the  clue  to  microcrack 
toughening.  Some  of  these  issues  have  been  addressed  by  Faber  [9]  and  we  will 
consider  such  issues  in  the  future. 

ACKNOWLEDGEMENT 

This  research  was  supported  by  the  Office  of  Naval  Research  through 
Contract  NR  064-N00014-81-K-0650  with  the  University  of  Illinois. 


-14- 


DISCUSSION 

Using  the  concept  of  the  mlcrocrack  density  e  first  introduced  by 
Budiansky  and  O'Connell  [1]  together  with  a  further  simplified  microcracking 
law  proposed  by  Fu  and  Evans  [10],  we  establish  the  basis  of  a  continuum 
mechanics  description  for  the  phenomenon  of  stress-induced  microcracking  in  a 
brittle  polycrystalline  composite.  With  appropriate  integration  of  the 
constitutive  relations  we  were  able  to  obtain  information  as  to  what  happens 
in  the  process  zone  at  the  vicinity  of  the  crack  tip  of  a  major  crack  under  a 
mode  I  loading  in  a  material  susceptible  to  microcracking.  For  appropriate 
choices  of  parameters,  the  transformation  zone,  as  in  the  case  of  small  scale 
yielding  in  a  ductile  material,  is  well  contained.  The  stress  and  strain 
fields  exhibiting  the  characteristics  shown  in  fig  (7)  and  (8)  are  well 
behaved.  When  the  first  critical  stress  <3q  and  the  saturation  value 
eM  are  specified,  the  second  critical  stress  Ojj  or  the  variable  X  appearing  in 
the  microcracking  law  control  the  rate  at  which  the  material  microcracks  and 
affects  the  stress  and  strain  fields  as  well  as  the  shape  and  size  of  the 
transformation  zone* 

At  this  point,  we  have  all  the  necessary  tools  for  a  full  scale 
investigation  of  the  phenomenon  of  microcracking.  Such  an  investigation 
should  address  the  question  of  the  possible  effects  on  material  toughening  due 
to  microcracking.  Further  work  is  necessary  to  elucidate  the  nature  of  the 
crack  propagation  criterion.  In  addition  studies  done  by  McMeeking  and  Evans 
[6]  and  Budiansky  et  al.  [7],  suggest  that  the  toughness  of  certain  ceramics 
can  be  substantially  enhanced  through  the  controlled  use  of  martensitic 
transformation.  Faber  [9]  found  that  microcracking  contributes  to  the 

toughness  of  Zr02  -  ceramics  and  Fu  and  Evans  [10]  found  that  stress  induced 
microcracking  enhances  the  toughness  of  brittle  polycrystalline  aggregates. 


the  microcrack  density  is  saturated  is  also  characterized  by  the  corresponding 
quantities  x8,  h8  and  08  fig.  (5). 

Fig.  (9)  to  fig  (9y)  show  how  the  above  quantities  depend  on  the  second 
critical  stress  Cjj,  for  fixed  values  of  oc  and  ■  X(<jjij-oc)  .  The  size  of 
the  fully  mlcrocracked  zone  is  inversely  proportional  to  the  value  of  the 
second  critical  stress  °M*  Similarly  high  saturation  values  of  the  microcrack 
density  yield  small  saturation  zone.  There  is  an  Increase  in  the  size  of  the 
above  zone  in  both  directions  as  X  increases  or  °M  decreases.  In  no  case  does 
the  inner  zone  becomes  Identical  with  the  whole  damaged  zone  because  as  we 
said  earlier  the  existences  of  a  zone  of  intermediate  microcracking  is  an 
essential  feature  of  the  microcracking  law  eq(3).  The  microcrack  zone 
boundary  depends  on  the  first  critical  stress  o^.  For  fixed  values  of  both 
the  first  critical  stress  oc  and  the  saturation  microcrack  density  the 
microcrack  boundary  tends  to  stretch  along  the  x  axis  while  reduction  of  h,,,  is 
observed,  as  X  Increases.  The  angles  0m  and  0g  corresponding  to  hm  and  hg 
drop  from  78°  for  small  values  of  X,  to  49°  as  X  becomes  infinite.  Fig  (10) 
shows  microcrack  zones  obtained  through  our  finite  element  analysis.  Fig  (11) 
and  (12)  show  the  profile  of  the  microcrack  density  ahead  of  the  crack  tip  for 
various  values  of  °M* 

Plane  stress  results 

Similar  results  to  those  of  the  plane  strain  case  were  obtained  for  the 
plane  stress  case.  The  transverse  constraint  of  the  plane  strain  case 
produces  large  transformation  zones,  the  overall  response  of  a  microcracking 
material  was  found  to  be  similar  for  both  plane  strain  and  plane  stress  cases. 
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Fig.  11  Microcrack  density  ahead  of  the  crack  //Jo 
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Fig.  12  Microcrack  density  ahead  of  the  crack  tip 


